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ABSTRACT 

Measurement of redshifted 21-cm emission from neutral hydrogen promises to 
be the most effective method for studying the reionisation history of hydrogen and, 
indirectly, the first galaxies. These studies will be limited not by raw sensitivity to 
the signal, but rather, by bright foreground radiation from Galactic and extragalac- 
tic radio sources and the Galactic continuum. In addition, leakage due to gain errors 
and non- ideal feeds conspire to further contaminate low-frequency radio obsevations. 
This leakage leads to a portion of the complex linear polarisation signal finding its 
way into Stokes /, and inhibits the detection of the non-polarised cosmological signal 
from the epoch of reionisation. In this work, we show that rotation measure synthesis 
can be used to recover the signature of cosmic hydrogen reionisation in the presence 
of contamination by polarised foregrounds. To achieve this, we apply the rotation 
measure synthesis technique to the Stokes / component of a synthetic data cube con- 
taining Galactic foreground emission, the effect of instrumental polarisation leakage, 
and redshifted 21-cm emission by neutral hydrogen from the epoch of reionisation. 
This produces an effective Stokes / Faraday dispersion function for each line of sight, 
from which instrumental polarisation leakage can be fitted and subtracted. Our re- 
sults show that it is possible to recover the signature of reionisation in its late stages 
(z w 7) by way of the 21-cm power spectrum, as well as through tomographic imaging 
of ionised cavities in the intergalactic medium. 

Key words: cosmology: observation — cosmology: theory — diffuse radiation - 
instrumentation: interferometers — radio continuum — techniques: polarimetric 



1 INTRODUCTION 

Statistical observations of the redshifted 21-cm emission 
from neutral hydrogen during the epoch of reionisation 
(EoR) promise to provide a wealth of information about 
the properties of neutral hydrogen at high redshift, as well 
as some of the fundamental astrophysics behind the reioni- 
sation process and the first luminous objects. While density 
perturbations in the matter distribution mediate fluctua- 
tions in the 21-cm signal both prior to and following reioni- 
sation, during the reionisation era the relation between the 
21-cm power spectrum and the underlying matter power 
spectrum is complex and, in its late stages, dominated by the 
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formation of large ionised "bubbles" ( Furlanetto et al.|2004 



McQuinn et al. 20061. These bubbles of ionised hydrogen 



imprint features on the 21-cm power spectrum that reflect 
the luminosity and clustering of ionising sources responsible 
for reionisation. The 21-cm power spectrum can be used as 
a statistical test to distinguish candidate reionisation mod- 
els, as well as to constrain the history and morphology of 
reionisation (Barkana 20091. 



Measuring this weak cosmic signal will be challenging 
however, owing to contamination from a large number of as- 
trophysical and non-astrophysical components. Indeed, as- 
trophysical foregrounds are typically brighter than the cos- 
mological 21-cm signal by 4-5 orders of magnitude (see, e.g., 
Pi Matteo et al.|2002|rOh fc Mack|2003| . Some of these fore- 



grounds also have bright polarised counterparts. As a result, 
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significant contamination will occur in the non-polarised 
(Stokes 7) component of observations due to non-ideal feed 
configuration and electronic errors in interferometers. This 
process is known as instrumental polarisation leakage. Var- 
ious experiments are planned to measure 21-cm emission 
from the neutral intergalactic medium (IGM), including the 
Low Frequency ArrajFj (LOFAR), the Murchison Widefield 
Arra>[^] (MWA) and the Precision Array to Probe Epoch of 
Reionizatiorj^] (PAPER). These instruments will be capable 
of operating at full-Stokes, meaning they will be capable of 
separating non-polarised signals from polarised signal. 

It has been suggested that foregrounds can be removed 
through continuum subtraction (e.g., |Gnedin fc Shaver|2004| 
Wang et al.|2006[ ), or using the differences in symmetry from 
power spectra analysis ( Morales & Hewitt 2004 , Zaldarriaga 



et al.|2004"|). This has been shown in recent work ( |Jelic et al 



2008 |Geil et al.|2008 1 to be an effective method for the pur- 



pose of detecting the EoR signal. On the other hand, elim- 
inating contamination of Stokes I due to instrumental po- 
larisation leakage has not been treated. This paper demon- 
strates a method that removes polarised contamination by 
utilising a procedure known as rotation measure (RM) syn- 
thesis. 

Rotation measure synthesis exploits the Fourier rela- 
tionship between the complex linear polarisation signal, 
P(A 2 ), and the Faraday dispersion function, F(4>), making 
it possible to identify components of polarised signal with 
specific Faraday depth, <j> (see, e.g., Burn|196 6 Brcnt jens fc| 
|de Bruyn||2005[ ). We modify this technique so as to make 
it applicable to the non-polarised Stokes I signal, calculat- 
ing an effective Stokes I Faraday dispersion function, Fi ((f)). 
The three-dimensional Stokes / data set is first re-sampled 
along the frequency axis into the natural A 2 -basis for Fara- 
day rotated emission. A one-dimensional Fourier transform 
is then applied along each line of sight. Any feature observed 
in F(4>) with a non-zero Faraday depth will, by the definition 
of our model, be polarised contamination, thereby making 
it possible to identify peaks of leaked polarised signal in 
Fi((f>). An iterative cleaning process, together with knowl- 
edge of the instrument's specific leakage behaviour, can be 
applied to Fi(<f>) in order to attempt to subtract this con- 
tamination. An inverse one-dimensional Fourier transform 
is then applied along each line of sight, followed by a coordi- 
nate transformation back to frequency space. Applying this 
procedure to all lines of sight in a field of view produces a 
three-dimensional Faraday dispersion cube; hence the use of 
the term Faraday tomography. 

This paper is organised as follows: Section [2] discusses 
the astrophysical sources we model in this work. Section 
discusses the basic principles of radio polarimetry. Section 
discusses the astrophysical and instrumental effects we simu- 
late in this work. The generation of synthetic Stokes / obser- 
vation data cubes is outlined in Section[5] Section [6] explains 
the continuum foreground removal strategy and introduces 
the leaked polarised foreground component subtraction algo- 
rithm we have developed. This section also includes a sum- 
mary of the assumptions used throughout. In Section [7] we 
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discuss the use of spherically averaged power spectra as a key 
statistical measure of cosmic reionisation, and our method 
of computing statistical instrumental sensitivity. Section [8] 
presents the results of our simulations and their analyses. 
Section [9] is a summary of this paper. 

Throughout this paper we adopt a set of cosmolog- 
ical parameters consistent with the Wilkinson Microwave 
Anisotropy Probe ( WMAP) ( |Komatsu et al.|2009| for a flat 
ACDM universe: Sl m = 0.27 (dark matter and baryons); 
Qa = 0.73 (cosmological constant); fib = 0.046 (baryons); 
h — 0.7 (Hubble parameter); n s = 1 (primordial spectrum 
index); erg = 0.8 (primordial spectrum normalisation). All 
distances are in comoving units unless stated otherwise. 



2 ASTROPHYSICAL SOURCES 

In this work, we model two astrophysical sources: the non- 
polarised 21-cm emission from the epoch of reionisation and 
diffuse Galactic synchrotron emission (DGSE). We consider 
the total intensity of DGSE which appears in the non- 
polarised Stokes I parameter, as well as the its polarised 
component. 



2.1 Epoch of reionisation signal 

Reionisation starts with isolated regions of ionised hydro- 
gen (Hn) forming around galaxies and quasars, which later 
grow and merge to surround clusters of galaxies. The reioni- 
sation process is complete when these regions overlap and fill 
the volume between galaxies. Stars are expected to be the 
main ionising sources during the EoR (see, e.g. , |Srbinovsky| 
& Wyithe||2007l iFaucher-Giguere et al. 1120081 |Loeb||2009h, 



although other sources such as quasars have also been sug- 
gested. To which level each of these contribute is still an open 
question. The stellar contribution begins with the formation 
of massive, metal-free stars that form in the shallow poten- 
tial wells of the first collapsed dark matter halos (z ~ 15), 
and continues by way of a normal population of less mas- 
sive stars that form from the metal-enriched gas in larger 
dark matter halos present at lower redshift (z ~ 6). We 
use a semi- numerical model for the reionisation of a three- 
dimensional volume of the IGM by galaxies. This model is 
described very briefly below and we direct the reader to |Geil| 
fc Wyithe| ( |2008| > for further details. 

Briefly, we begin by simulating the linear matter over- 
density field S(x, z) = pui(x, z)/p m — 1 inside a periodic, co- 
moving, cubic region of volume V = L 3 , by calculating the 
density contrast in Fourier space 5(k, z) corresponding to 
a ACDM power spectrum (Eisenstein & Hu 19991 extrapo- 
lated to a specified redshift z. The semi-analytic model used 
to compute the relation between the local dark matter over- 
density and the ionisation state of the IGM is based on the 



model described by Wyithe fc Loeb (20071 and Wyithe fc 



Morales (20071. The sum of astrophysical effects included 



in our model leads to the growth of H n regions via a phase 
of percolation during which individual Hu regions overlap 
around clustered sources in overdense regions of the uni- 
verse. Where necessary, we consider a model in which the 
mean IGM is fully reionised by redshift z ov = 6 ( Fan et al. 
2006}|Gnedin fc Fan|2006| |White et aTp553 i. 
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Figure 1. Three-dimensional model epoch of reionisation signal at z = 7 showing fully ionised regions with no 21-cm emission (blue, 
mK) and partially ionised, 21-cm signal emitting regions (red, 11 mK). Each box has a comoving side length of approximately 524 Mpc. 
Left: Non-instrumentally smoothed. Centre: Instrumentally smoothed in the sky-plane using a Guassian tapered synthesised beam with 
a 5 arc minute full-width half-maximum. Right: Instrumentally smoothed, central frequency channel slice. 



The corresponding average 21-cm brightness temper- 
ature contrast of a gas parcel with the cosmic microwave 
background (CMB) ( Madau et al.|1997[ ) is given by 



ST b 



26 x Hl (l + S b ) 1 



T 



0.022 



1/2 



> l^ 15 ,^] mK. 

Qmh 2 10 



(1) 



Here, the spin temperature, T s , CMB brightness tempera- 
ture, T^, baryonic overdensity, 8b, (which we assume to be 
equal to the dark matter overdensity) and neutral fraction, 
xh i , are calculated for a location x at redshift z. We as- 
sume T a T 7 during the epoch of reionisation (Ciardi & 
Madau 2003 Furlanetto et al. 2006 ) and ignore the enhance- 



ment of brightness temperature fluctuations due to pecu- 
liar velocities in overdense regions ( |Bharadwaj fc Ali||2005 



Barkana fc Loeb|2005 I. Peculiar velocities were included in 
the semi- numerical model of Mesinger & Furlanetto (20071, 



who found their effect to be small on scales ~ 10 Mpc. Fig- 
ure [l] shows an example of a synthetic three-dimensional 
epoch of reionisation signal data cube at z = 7 with and 
without the effect of limited instrumental angular resolution 
(see Section 4.2 for a description of the synthesised beam 
resolution). 

For the purpose of demonstrating foreground removal, 
we do not consider the evolution of the IGM in the line-of- 
sight direction within the observed bandpass. Therefore, the 
realisation of the density and ionisation state of the IGM in 
each frequency channel is considered to be at the same stage 
of cosmic evolution. This is a valid approximation when the 
range in redshift corresponding to the simulated bandwidth 
is small compared with the central redshift of the bandwidth 
Az/z <C 1 for redshifts greater than the redshift of overlap 
z > z ov . We also assume that the reionisation signal is non- 
polarised and is therefore present in Stokes I only. 



2.2 Diffuse Galactic synchrotron emission 

Diffuse Galactic synchrotron emission originates from rela- 
tivistic electrons in the interstellar medium interacting with 
the Galactic magnetic field. Although we are primarily con- 
cerned with the polarised component in this paper, we dis- 
cuss our model for total intensity of synchrotron emission 
here because of the proportionality between the two, as well 
as to demonstrate the combined effect of continuum and 
polarised-component foreground removal. For more detailed 
non-polarised and polarised foreground models, including 
synchrotron emission from discrete sources such as super- 
nova remnants and free-free emission from diffuse ionised 
gas, see 
(20101 



Jelic et al. (20081; Bowman et al. (20091; Jelic et al 



2.2.1 Total intensity of diffuse Galactic synchrotron 



emission 



Foreground contamination and its removal could have sig- 
nificant consequences for the detectability of the weak, red- 
shifted 21-cm signal. Indeed, the total intensity of fore- 
grounds will be brighter than the cosmological 21-cm sig- 
nal by 4-5 orders of magnitude. The three main sources 
of foreground contamination of the 21-cm signal are DGSE 
(which comprises ~ 70 per cent near 150 MHz), extragalactic 
point sources (~ 27 per cent) and Galactic bremsstrahlung 
(~ 1 per cent) (Shaver et al. 19991. The frequency depen- 



dence of each of these foregrounds can be approximated by 



a power law with a running spectral index (Shaver et al. 
1999| |Tegmark et al.||2000[ ). While the sum of power laws is 
not in general a power law, over a relatively narrow fre- 
quency range, B, (such as that considered in this paper 
where B/v <C 1), a Taylor expansion around a power law 
can be used to describe the spectral shape. We therefore also 
approximate the sum of foregrounds as a power law with a 
running spectral index, and specialise to the case of Galactic 
synchrotron emission, which dominates the foregrounds. For 
the purpose of demonstrating the removal of polarised fore- 
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v = 177.5 MHz 




6 X [arc minute] v [MHz] 

Figure 2. Left: A single frequency slice at v = 177.5 MHz showing the angular structure of the diffuse Galactic synchrotron foreground 
emission, for which (T^,a T ) RJ (183, 16) K. Right: DGSE as a function of frequency for a number of random lines of sight. The dashed 
line shows the value of T^(v) for ST^ = 0. 



ground leakage, we assume that the brightest point sources 
have already been removed. 

The total intensity of Galactic synchrotron emission 
varies as a function of both sky position and frequency. 
We model the frequency and angular dependence of Galac- 
tic synchrotron foreground emission (following |Shaver et al.| 
1999| |Tegmark et al.||2000| |Wang et al"]|2006[ ) by first con- 
structing a realisation of the angular fluctuations in the 
foreground, <5T G (#), at a particular frequency within the 
band of interest using a two-dimensional Gaussian random 
field. We then add a mean sky brightness, T G , to our two- 
dimensional realisation of brightness temperature fluctua- 
tions. This foreground plane is then extended into three 
dimensions by extrapolation, using the running power law 
form along each line of sight, 



T G (0>) 



[f£ + ST£ (9)\ 



1A, 



(2) 



We have modelled the sky-plane variation of a syn in this 
paper by drawing a value for Aa syn from a Gaussian proba- 
bility distribution with a standard deviation of 0.1, for every 
line of sight. Note that over a data cube at the frequencies 
of interest, T G ~ 200 K. The strategy for removing contin- 
uum foreground contamination is discussed in Section |6.1| 
Figure [2] includes a single frequency slice at v = 177.5 MHz 
(left), showing the angular structure of the total intensity of 
diffuse Galactic synchrotron foreground emission, for which 
(T b G ,o- T ) = (183, 16) K. Also shown (right) is the DGSE as 
a function of frequency for a number of lines of sight. The 
dashed line shows the value of T^(v) for 8T^ — 0. 

Although the mean brightness temperature of the 
DGSE is ~ 200 K, the interferometric response to this emis- 
sion will only be the fluctuation from the mean brightness 



temperature within the primary bearrj^] This is because the 
signal recieved by any single antenna does not contribute 
to the synthesised signal, i.e., the instrument as a whole re- 
sponds to non-zero baseline spacing only. Such fluctuations 
are generally much smaller than the mean over the primary 
beam. We retain the primary beam's mean brightness tem- 
perature in our DGSE model with no consequence, since any 
constant signal along a line of sight is removed during the 
continuum foreground subtraction stage. 



2.2.2 Polarised DGSE 

The polarisation state of a stream of radiation may be spec- 
ified by the four Stokes parameters /, Q, U and V. Linearly 
polarised radiation can be specified using Q and U only and 
is often represented by the complex linear polarisation P, 
given by 



P = Q + iU =plexp(2ix), 



(3) 



where Q = 5ft [P] and U — Q[P], p is the intrinsic degree 
of polarisation and \ ls the position angle (or polarisation 
angle) of the linearly polarised radiation. These quantities 
are related by 



P i = VQ 2 + u 2 , x 



■ arctan 



U 



(4) 



We model the polarised component of the DGSE as be- 
ing proportional to its total intensity, T^(v), such that 



Q G (u) = pT b G (i/)cos(2x), 
U G (u) = pT b G Hsm(2 X ), 



(5) 
(6) 



4 The primary beam is the field of view of a single antenna tile. 
MWA antennas have fields of view of ~ 240A 2 square degrees (to 
the first null), where A is in the range 1 < A < 4 m. 
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where T^{u) is the total Galactic synchrotron bright- 
ness temperature, x is the polarisation angle (line-of-sight- 
dependent) and the polarisation fraction, p, is dependent on 
the brightness temperature spectral index, a syn , through 



P = 



3a s 



0.7, 



(7) 



( Le Roux|196l l 



|Bernardi et al. ( 2009 2010 1 give constraints on the level 



of polarised DGSE using Westerbork Synthesis Radio Tele- 
scope data at 150 MHz and arcmin scales. They found that 
the rms value for polarised foreground is « 7.2 K on 4 arc 
minute scales. Using a polarisation fraction of 0.7, the fluc- 
tuations in our polarised foreground at the same frequency 
with a 5 arc minute full width half maximum beam are of 
the same order of magnitude (0.7 x 25 K = 18 K). 



3 BASICS OF RADIO POLARIMETRY 
3.1 Faraday rotation 

The rotation of the plane of polarisation of radiation (spec- 
ified by x) propagating through a magneto-ionic medium 
results from the birefringence of the plasma. This effect, 
known as Faraday rotation, has a nonlinear dependence on 
frequency such that 



X = Xo + 



(8) 



where <j) is the Faraday depth and xo is the non-rotated 
polarisation angle. Faraday depth is defined by 



8.1 x 10 a 



,B ■ dr rad m 



(9) 



where n e is the electron density in cm" , B is the mag- 
netic field in Gauss, L is the distance between the point 
of emission and the observer in parsec and dr is an in- 
finitesimal path vector in parsec (away from the point of 
emission). A positive Faraday depth indicates an integrated 
magnetic field pointing toward the observer, however, the 
Faraday depth may be negative due to the integrated orien- 
tation of the magnetic field pointing away from the observer 
along the line of sight. 

Along any particular line of sight, there may be any 
number of Faraday screens at different Faraday depths by 
which polarised radiation may be Faraday rotated having 
been emitted with a unique initial polarisation angle. We 
include only one screen per line of sight in our simulations 
in order to show the principle behind the instrumental po- 
larisation leakage cleaning process. 

There are two planned target EoR fields for the MWA, 
centred on (a, 8) = (4 h , -30°) and (7 h 40 m ,10°), which cor- 
respond to the Galactic coordinates (I, 6) « (228°,— 49°) 
and (209°, +15°) respectively. The interpolated RM values 
(and thus Faraday depth) of both of these field s including 
extragalactic sources are 20 <C 4> ^ 50 radm -2 (Johnston- 
Hollitt et al.|2004 1. It is reasonable to assume that tiie mag- 
nitude of Faraday depth is either less than or not much 
larger than extragalactic RM values in the same direction. 
RM ranges i n other directions have been giv en by|Haverkorn 



et al. 



(2010) 



and 



Bernardi et al. 



(2003|) (-17 <; 4> <; 10 radm 

<; 10 radm" 2 ). In our synthetic data cube, we 



have chosen to include a single Galactic Faraday screen with 
1 <C <f> <J 5 radm -2 in every line of sight. This range of 
shallow Faraday depth has been chosen despite the larger 
observed values since we have found contamination by po- 
larised signal passing through screens at small Faraday 
depth more troublesome to clean. Therefore, our results rep- 
resent a worst case scenario. We model the angular distri- 
bution of Faraday depth of the Galactic screen by a small- 
gradient-plane over the sky-plane (such that, for generality, 
V(j>(x, y) oc x + y). This is done to approximate the small, 
slowly varying Faraday depth over the proposed EoR fields 
of the MWA. 



3.2 Faraday dispersion function 

In dealing with the RM synthesis of a distribution of 
synchrotron-emitting and Faraday rotating regions along a 
line of sight, we follow the work of Brent j ens fc de Bruyn| 
(2005) and direct readers to their paper for a more thor- 



ough explanation, as well as Schnitzeler et al. ( 2009 1 for an 
application of the technique. 



3.2.1 Continuous case 

The Faraday dispersion function, F(4>), can be defined 
through 



P(A 2 ) = / F(rj))exp(2i(f)\ 2 )d(j>, 



(10) 



where the frequency dependence of P has been given in 
terms of A 2 since A 2 and <f> form a Fourier-variable pair. In 
practice, the sampling of A 2 will be limited by the finite fre- 
quency bandpass of the correlator. Therefore, we introduce 
a weighting or sampling function, W^A 2 ), such that 



Po ba (A 2 ) = W{\ 2 ) / F(0)exp(2i<M 2 ) 



where, for example, 
W(\ 2 ) = 



if channel is recorded 
if otherwise. 



(11) 



(12) 



Other forms of the sampling function may include non- 
uniform weighting to account for non-ideal response for cer- 
tain channels (at the shoulders of the bandpass, for exam- 
ple). Equation ( 11 1 can be Fourier inverted to find F(cf>) and 
convolved with the inverse Fourier transform of the sam- 
pling function. This gives the Faraday dispersion function 
that would be obtained given the limited sampling of P, 

F ohs ((t>) = F(<j}) * R(4>) 

= K j Po bs (A 2 )exp(-2^A 2 )dA 2 , (13) 

where the rotation measure spread function, R((f)), is defined 
by 



R(<j})=K j W(A 2 )exp(-2i0A 2 )dA 2 , 
and the appropriate normalisation factor is 
K = I W(X 2 ) dX 2 



(14) 



(15) 
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Figure 3. Demonstration of requirements for back rotation of 
polarisation angle showing the magnitude {solid black), real part 
{dashed red) and imaginary part {dashed blue) of R{<f>) for a 
32 MHz bandpass of uniform response centred on 177.5 MHz. 
Bottom: Polarisation angle derotated using A 2 , = 0. Top: Polari- 
sation angle derotated using the weighted mean value of A 2 . 



The rotation measure spread function is a complex- 
valued function whose magnitude is peaked and symmetric 
about 4> = 0. The lower panel of Figure [3] shows the mag- 
nitude, real part and imaginary part of R(<f>) for a 32 MHz 
bandpass of uniform response centred on 177.5 MHz. It can 
be seen that both the real and imaginary parts oscillate 
rapidly throughout the central response peak. This makes 
it difficult to determine x> which given the form of Equa- 



tions ( 13 1 and ( 14 1 is the polarisation angle derotated back 
to A 2 = 0. It is possible, instead, to shift the A 2 variable 
by some amount, affecting the oscillitory behaviour of the 
complex components of R(4>). In general, the optimal shift 
is given by the weighted mean value of A 2 , 



Aq = / W(A 2 )A 2 dA 2 



W(\ 2 )d\ 2 , 



(16) 



It is simple to show, as can be seen in the top panel of Fig- 
ure [3] that this forces dQ[R{(p)]/dtf>\ 4>=o = 0, thereby making 
the evaluation of the polarisation angle using Equation Q 



more accurate. Hence, using Equation (jsj, we can now write 
the derotated polarisation angle as 



Xo = x(-V>) - <Mo- 



(17) 



Dropping the 'obs' subscript hereafter and making the 
shift in A 2 , we have 



P(A 2 ) = W(A 2 ) J F(4>)exp [2i(f>{\ 2 - A 2 ,)] d$ 
F{4>) = K J P(A 2 )exp [~2icj>{\ 2 - A 2 ,)] d\ 2 , 
R{4>) = K [ W{\ 2 ) exp [-2#(A 2 - A 2 ,)] dA 2 



(18) 
(19) 
(20) 



where K is defined as before in Equation ( 15 1 



3.2.2 Discrete case 

Since a correlator samples signals discretely, we must use 
discrete rather than continuous Fourier transformations. To 
do this, consider a correlator which provides iV c h frequency 
channels over a bandpass of width B centred on vq, such 
that each channel has a width Ais = B/N c t. The correlator 
gives signal measurements for equally spaced channels in 
frequency, however, the parameter we wish to work with is 
A 2 . We make the following approximations for the channel 
centres and widths. The central A 2 of channel j is given by 



A 



/■ 


A 2 di^/ 


7 






c 2 




( Ai/ 


-I 





dv 



+ o 











(21) 



where Vj is the central frequency of the channel. The width 
of channel j is given by 



AA 2 



1 + 



1 



+ o 




[(f) 4 ] 







(22) 



If 0AA < 1 for all channels, Equations ( 15 1, ( 19 1 and ( 20 1 



can be approximated by the following Fourier sums: 

JVch-i 

F(0) w K Y, p i ex P [ — 2i0(A| — A 2 ,)] AA 2 ; 

3=0 
N ch -1 

R{4>) w K Y w i ex P [~2i<j!>(A 2 — Aq)] AA 2 ; 

3 = 



(23) 
(24) 
(25) 



where A 2 is the central A 2 of channel j as approximated by 



Equation (21 1 



In the lag-domain, the correlator we consider corre- 
sponds to iV c h lag channels covering a total time At = 1/Af, 
with lag spacing St = 1/B. By choosing to perform the 
frequency-lag conversion with respect to the central fre- 
quency, the recorded lag spectrum will be centred on t = 0. 
This enables the recovery of Faraday depths of either sign 
in the Faraday dispersion function. 
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3.3 Equivalent Faraday dispersion function for 
Stokes I 

Unlike the complex linear polarisation, P, the non-polarised 
Stokes I signal is a real-valued function. Despite this, it is 
possible to calculate an effective Stokes / Faraday dispersion 
function for this signal by assigning its imaginary part to 
zero in all channels. Then, in analogy with Equation (23 1, 



Fi(<j>) » K 



3=0 



Ij exp [— 2i(f>(\ 



Ao) 



AX1 



(26) 



where Ij is the Stokes I signal measurement for channel j. 
Since the Stokes I response is a real function, by the prop- 
erty of the Fourier transform, the effective Stokes I Fara- 
day dispersion function is therefore Hermitian, i.e., Ff(<j)) = 
Fi(—(j>). Because of this redundancy, we may plot Fi(<f>) for 
<j> ^ only without any loss of information. 
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Figure 4. Bandwidth depolarisation for the MWA [Ai> = 40 kHz 
centred on v c = V2\/(l + z)\. 



4 ASTROPHYSICAL AND INSTRUMENTAL 
EFFECTS 

In this section, we discuss some of the astrophysical and in- 
strumental effects relevant to radio polarimetry and radio 
interferometry in general, and introduce the effect of instru- 
mental polarisation leakage. For convenience, we consider 
the propagation of the Stokes vector visibility throughout 
the system in terms of the brightness temperature. 

Throughout, we consider a radio interferometer with 
specifications comparable with the MWA. The MWA will 
consist of w 500 tiles, each with 16 cross-dipoles with 1.07 m 
spacing and an effective area A E w 16(A 2 /4) (A c is capped 
for A > 2.1 m). For simplicity, we model the beam-formed 
response of each tile as that of a single cross-dipole. We as- 
sume a continuous antenna density profile, p a oc r -2 , within 
a 750 m radius with a const ant c ore density of approximately 
one tile per 36 m 2 . Section Al of the Appendix gives more 



detail on the fiducial array configuration considered. 



Stokes Q and U spectra and then summing the contribu- 
tions over each channel to produce the final spectra with 
the desired number of channels. 

Note that the detectable range and sampling resolution 
of Faraday depth can be found in the following way. The 
sinc-f actor effectively enveloping P c ^ vanishs when a — tt/2. 
Solving for <f> when P c h = gives 



4A 2 Ai/ 



(29) 



where A 2 , is the corresponding central value of A 2 given by 



Equation (16 1. From this, the Faraday depth sampling in- 
terval is 



A(f): 



20n 



Nr.: 



2A 2 P' 



(30) 



Note that the central lag channel (if N c h is odd, otherwise 
the central two channels) contains all signals with \(f>\ < A</>. 



4.1 Bandwidth depolarisation 

The integrated complex polarisation over each frequency 
channel will suffer from depolarisation due to the decon- 
structive summation of linear polarisation vectors within the 
small but finite range of frequencies in each channel. The in- 
tegrated complex polarisation for a channel with frequency 
range v\ ^ v ^ v-x is given by 



P(u) du. 



dv. 



(27) 



v\ centred on v c 



Given a channel of width Av = U2 
(y\ + V2)/2, such that Av/v c <C 1, 

. . sin a 2i(t>\ 2 rr , s 

P c hM=P e v <=, (28) 

a 

where a = 2(f)\^.Av/v c , and it is assumed p is independent 
of frequency. Figure [1] shows the bandwidth depolarisation 
for the MWA centred on v = 143 MHz and 204 MHz (cor- 
responding to the redshifted 21-cm emission frequency for 
z — 9 and 6 respectively). We include the effects of band- 
width depolarisation by generating over-sampled synthetic 



4.2 The synthesised beam and beam 
depolarisation 

Due to the finite length of the longest baseline, radio in- 
terferometers have a correspondingly finite angular reso- 
lution. This can be quantified in terms of the synthesised 
beamwidth, given by 



A6 h = 



0.705A 
D ' 



(31) 



where A is the wavelength of observing and D is the maxi- 
mum baseline lengtrj^] When calculating the actual angular 



5 Note that the beamwidth referred to by many authors is A9 = 
1.22A/D ~ A/D, which is the highest angular resolution meet- 
ing the Rayleigh criterion. This corresponds to the central peak 
to first null separation in the resulting beam-response pattern of 
a circular aperture. Other conventions include: full-beamwidth 
between first nulls, A^f^vfn = 2.44A/D; full-beamwidth at 
half-power, A6>fwhp = 1.02A/D; and full-beamwidth at half- 
maximum, A6pwHM = 0.705A/D. We have adopted the last of 
these conventions in this paper. 
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resolution achievable in the field of observation, the pro- 
jection of the baselines onto the sky-plane must be taken 
into consideration. Throughout this work we assume obser- 
vations of a zenith-field only. 

In addition to the limited angular resolution imposed 
by the maximum baseline length of the interferometer, sky- 
plane images can be smoothed by applying a weighted filter 
to the two-dimensional Fourier transform of each image for 
every frequency channel. This can be used to take weight 
from data that have been poorly sampled due to insuffi- 
cient baseline coverage for example, at the expense of angu- 
lar resolution. The particular type of filter used depends on 
the final observation product required. For all 21-cm sim- 
ulations in this work we have applied a Gaussian filter in 
Fourier-space which corresponds in image-space to a Gaus- 
sian smoothing function of width of 5 arc minutes for all 
frequency channels. 

In addition to bandwidth depolarisation, a second cause 
of depolarisation is beam averaging. Unresolved polarised 
sources within a synthesised beam will add vectorally and 
may therefore superimpose destructively. Similarly, a spatial 
gradient in Faraday depth over the synthesised beam will 
produce a variation in the angle of polarisation as a func- 
tion of position. If this variation is significant when averaged 
across a single synthesised beam, the polarised fraction can 
be reduced. For simplicity, we have not modelled these ef- 
fects in this work. 



Jones operators ( Jones||l941 1, J, transforming as 

= Ji (33) 

where the prime indicates the vector after transformation. 

The combined effect of these transformations can be 
mathematically described by a transformation made up of 
the product of any number of Jones matrices. For simplicity, 
we only consider Faraday rotation, complex gain and polar- 
isation leakage in this work, assuming that all other effects 
have been accounted for, by calibration for example. There- 
fore, the combined Jones matrix for interferometer arm 'i' 
is given by 



Ji — Gi Di Fi 



(34) 



where (in order of operation) Fi is the Faraday rotation 
operator, Di is the operator that accounts for polarisation 
leakage through deviations of the actual feed from the ideal 
and Gi is the complex receiver gain operator. These opera- 
tors are given by 



F, = 



Di 



cos Axi — sin A\i 
sin Axi cos Axi 



1 



l 



(35) 



(36) 



4.3 Representation 



In this section, we closely follow t he notation of |Hamaker| 
et al.| ( |1996[ ), |Sault et al.] ( |1996| and |Thompson et a.l.| ( |2001| ) . 

The electric field of the most general homogeneous plane 
wave propagating in the direction k, recieved at the feed po- 
sition of antenna V can be given with respect to the Carte- 
sian coordinates x, y and z by 



e*(t) = 



Ej.. x 

Ei,y 



exp[ 



Jt], 



(32) 



where the amplitudes Ei tX and Ei tV are complex numbers, 
allowing for a possible phase difference. If these components 
have the same phase the wave is said to be linearly polarised; 
if different it is said to be elliptically polarised, with the spe- 
cial case of a n/2 phase difference giving a circularly po- 
larised wave. 

The plane waves recieved by the antennas may have 
been modified by a number of spurious effects that occur 
along the path between the source and the feeds. The result- 
ing electrical signal is also modified by a number of inten- 
tional and other spurious effects that occur along the path 
between the feeds and the correlator. These effects include: 
physical effects such as the Faraday rotation by magneto- 
ionic media between the signal source and the feeds; geo- 
metrical effects such as the rotation of the feeds with re- 
spect to the polarisation coordinates of the signal, including 
the projection of the dipoles onto the observed sky plane; 
instrumental effects such as the conversion from an electric 
field to a voltage in the feed electronics; signal amplifica- 
tion; and the effects of polarisation leakage or variation of 
gain with frequency. These effects can be described by lin- 
ear transformations and represented by what are known as 



Gi = 









9i,y 



(37) 



where Axi is the angle of Faraday rotation, di. xy is the spu- 
rious sensitivity of the x receptor to the y polarisation and 
vice versa, and Qi. x is the gain error coefficient for the x re- 
ceptor. Note that Fi is real-valued whereas Di and Gi are 
complex- valued . 

For a single interferometric baseline, the measured vis- 
ibilities can be described by a 2 x 2 coherency matri^ 
Cij, which is formed by the time-averaged outer product 
of the two-dimensional, complex vector input signals ei and 
ej from the two antennas making up the baseline, 



; ei (*)®et(i)}, 



(38) 



where the dagger symbol represents the adjoint operation 
of complex conjugation and transposition. It is also possible 
to describe the measured visibilities by a four-dimensional 
coherency column vector, 



(ei(t) ®ej(*)) 
/ Ei. x E* x \ 



TP* 



E- E* 
V E i,v E j,v ) 



(39) 
(40) 



6 Throughout this paper, vectors and operators in the coherencey 
or Stokes representation are denoted by column vectors (lower- 
case bold non-italic) and matrices (upper-case bold non-italic) 
respectively, in order to distinguish them from vectors with no 
specific basis (bold italic). 
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Including the effect of the transformations, we have 

c'ij = (e'i(t) ®e'?(t)) 

= (J.e^t)® [Jje, (*)]*) 

= (Jiei(t) ® J*e*(t)} 

= (Ji ® Ji)(ej(t) ®e*(t)} 



(41) 



where Jy = Ji ® J* and we have made use of the following 
identity for the outer product: 

(AB)®(CD) = (A«C)(B®D). (42) 

As the Stokes parameters are used almost universally 
in radio polarimetry, we now make the conversion from the 
coherency representation to the Stokes visibility representa- 
tion. A distinction should be made at this point: the Stokes 
parameters, I, Q, U and V, are used when describing the 
intensity or brightness of radiation, whereas the Stokes visi- 
bilities are used to represent the complex visibilities as mea- 
sured by an interferometer. The context in which these terms 
are used provides the means for disambiguation. The Stokes 
visibility vector as measured by a single baseline made up of 
antennas i and j, By = [hj Qij Z7y Vij] T , is related to the 
coherency vector via a simple coordinate transformation, S, 
through 



S Ci 



where 



/ 1 








1 


1 








-1 





1 


1 





V o 


—i 


i 






(43) 



(44) 



the inverse of which is 



/ 1 1 \ 

1 j 

1 -i 

\ 1 -1 / 



(45) 



Note that the complex vector input signals and Jones matri- 
ces are two-dimensional whereas the coherence vectors and 
Stokes visibility vectors are four-dimensional. 

Operators that act in coherency- or Stokes-space can be 
transformed between representations using the same trans- 
formation matrix, S, used for the coherency and Stokes vec- 
tors, through 

A (S) =SA (C) S~\ (46) 

where is the operator in coherency representation and 
A' s ' is the operator in Stokes representation (operators in 
the Stokes representation that are used to reproduce the ef- 
fect of one or more non-depolarising, transforming elements 
are called Mueller matrices). Therefore, we can write the full 
Mueller matrix for the single baseline system as 



M, 



S(J 2 



.J*) 



(47) 



4.4 Thin Faraday screens 



Assuming the signal through both interferometic arms expe- 
rience the same degree of Faraday rotation Ax (a valid as- 
sumption for baselines shorter than the scale of fluctuations 



in the ionosphere), applying the coherency-Stokes visibility 
transformation on the Faraday rotation operators of a single 
baseline gives 



Mi 



/ 1 




V o 





cos(2Ax) 
sin(2A X ) 




\ 

-sin(2A X ) 

cos(2Ax) 

1/ 



(48) 



This is equivalent to a rotation of the P = [Q,U] T vector 
through an angle of 2Ax radian. So we see that Faraday ro- 
tation in the reduced Stokes representation (Q and U only) 
is the same as that in the complex vector representation 
(where F acts). This is reconcilable with Equation (J3J) since 
P is rotated by 2Ax = 2(j>X 2 from an initial angle 2%o re- 
sulting in a final angle of 2\ = 2%o + 2Ax- 

The efffect of thin Faraday screens on incident syn- 
chrotron emission can be modelled by 



k=i 

UF(\f) = ^pfc(Aj)sin(2xjfc). 



(49) 
(50) 



where Nf is the number of thin Faraday screens and \jk = 
(Xo)fc + 4>k{Xj — Ao), (xo)fc being the intrinsic polarisation 
angle at the fcth screen. 



4.5 Instrumental polarisation leakage 

In practice, feeds are non-ideal and this leads to leakage be- 
tween all Stokes parameters, I, Q, U and V. The treatment 
of instrumental polarisation leakage is similar to that as for 
Faraday rotation, but complicated by the fact that in gen- 
eral the feed errors corrupt all four Stokes parameters and 
cannot be assumed to be the same in all antennas. 

Assuming equal weights for all Stokes visibilities gives 
the following Mueller operator for the full interferometric 
system, 



M 



\: .\ i;, 

1 N N 

iV W^i)^5 s(J< ® j;)s " 1 ' 



(51) 



where both the baseline and its conjugate baseline are 
summed over in order to give a real- valued result. Using 
Equation ( |47| l, the Mueller matrix describing leakage for a 
single baseline with parallel linear feeds (i.e., cross-dipole 
feeds in antenna i are approximately parallel to those in 
antenna j) to first order in leakage and gain error is given 
by 



/ 



M, 



7++ 
7+- 
S+- 
-i5 



+ 



7+- 
7++ 
S++ 
iS— 



5+- 
S++ 

7++ 
-i7_ _ 



-i8- 
-i5- 



ij 

7++ / 



(52) 
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where I4 is the 4x4 identity matrix, 

7++ = (Aflv + Ag iiV ) + (&gix + A 9j,y), (53) 

7+_ = (A gi , x - Ag i>y ) + (Agl x - A ff * H ), (54) 

7— = ( A 9i,x ~ Ag it y) - (Agl x - Agl y ), (55) 

S++ = {di, xy + di,y X ) + (dj >xy + d* tVX ), (56) 

<5+— — (d-i,xy — di : y X ) + (dj tXV — dj yx ), (57) 

<5-+ = + <fc,vaO - ( rf jUy + (58) 

6 — = (di. X y di r y X ) {dj jXy dj yx ), (59) 

and uniform, unit ideal gain^] have been used for conve- 
nience, such that 



g%,p = 1 + Agi tP . 



(60) 



The resulting signal (again, for a single baseline only) includ- 
ing leakage and gain errors (primed) that originates from a 
general input signal (non-primed) is 



± I 1 



(1 + o7++) j *j + o7+-<3« + n8+-Uij - -iS-+Vij, 



(61) 



Q'ii = (1 + ^7++)Qij + \l+-^J + 2 S ++ U v ~ 5 < 



(62) 

(1 + ^7++)^ + o*+- / « _ o^+Qw + o*7 — 



(63) 

Figure [5] shows simulated examples of the Faraday disper- 
sion function, F(<j>), for a single baseline, simply illustrat- 
ing the effect of leakage between Stokes Q and U, noise 
and DC offsets. In these examples, gain errors and leakage 
from Stokes / have been ignored. The amount of leakage 
from Stokes Q into Stokes U and Stokes U into Stokes Q 
is parameterised by Lqu and Luq respectively (see Ap- 



pendix A3 1. The error coffecients in these example have been 



exaggerated for demonstration. The upper plot shows the 
three Faraday components only, at depths cj> = —100, 250 
and 300 radm -2 . The sidelobe structures that are visibly 
associated with each peak originate from the finite band- 
width of the correlator and are scaled copies of the rotation 
measure spread function, R(<j>), introduced in Section 3.2 



The relative magnitude of each peak depends on the intrin- 
sic brightness temperature of the linearly polarised radia- 
tion that originates from behind and between the Faraday 
screens. The absolute magnitude of the peaks are depen- 
dent on the size (and sign) of the leakage and gain errors 
for that particular baseline. The lower plot shows the same 
Faraday components (labelled 'F'), however, noise (labelled 
'N'), leakage (labelled 'L') and DC offsets (labelled 'DC') 
have been included. Power in the Faraday dispersion func- 
tion due to noise appears in all Faraday depth channels. Note 
that unlike the effective Stokes I Faraday dispersion func- 
tion, F((f>) is not Hermitian and therefore |F(— </>) 7^ F(<^)|. 

When considering all baselines, we simulate both gain 
and leakage errors by sampling a zero-mean, Gaussian dis- 
tribution. The standard variation of these distributions have 

7 For non-uniform and/or non-unit gains the relevant gain fac- 
tors follow through the algebra, skewing the Mueller matrix. This 
can be incorporated through the original choice of basis, leaving 
Equations l|52||-(|59| unchanged. 
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Figure 5. Examples of Faraday dispersion functions for a single 
baseline. Top: Faraday components only. Bottom: The same Fara- 
day components (F) with noise (N), leakage (L) and DC offsets 
(DC) added. 



been chosen to be 10 per cent of the nominal, unit gain for 
each polarisation and for each antenna. 



4.6 Instrumental calibration 

We assume that a real-time calibrating procedure will pro- 
duce a full polarisation calibration solution, including fits 
from polarised calibrators and constraints of most of the 
degrees of freedom in the primary beam models using non- 



polarised calibrators (Mitchell et al. 2008 Lonsdale et al 



2009). However, these calibration solutions will be produced 
from real-time analysis, in image-space and on short time 
scales (~ 10 s). It can be assumed these solutions will only 
be stable over some minutes and not the total integration 
time of hundreds of hours spread over months/years. In ad- 
dition, calibration solutions are to be applied on a baseline 
by baseline basis in order to obtain good polarisation purity 
and to prevent calibration errors from introducing artifacts 
into the snapshots. The instrument polarisation we address 
in this work, however, are the leakages and errors that do 
make their way into the snapshot observations and there- 
fore the fully integrated EoR data sets, that may still be 
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significantly larger than the cosmic signal. Therefore, using 
a single calibration solution is valid only for these short snap- 
shots, where the S/N is very low, thus hiding the polarised 
foreground contamination in a forest of peaks (in all Stokes 
parameters) . We overcome this by using a single calibration 
solution which we assume to represent some linear combi- 
nation of known calibration solutions from all snapshots. 

4.7 Instrumental noise 

Thermal noise is present in every visibility measurement for 
every frequency channel and will be instrument dependent. 
Section IA2.1I discusses how we model thermal noise. Since 
we are applying cleaning procedures to fully integrated data 
(as opposed to short snapshots), the rms of the instrumen- 
tal noise is generally much smaller than the contamination 
by polarised foregrounds considered in this work. As a re- 
sult, the fitting of peaks in Faraday dispersion functions is 
only slightly affected. We have found this to be the case 
when instrumental noise has been included in cleaning sim- 
ulations. For this reason, we do not include thermal noise in 
our simulations, but rather, provide an instrumental sensi- 
tivity curve when considering the statistical detectability of 
the EoR signal (see Section [71. 



5 SYNTHETIC DATA CUBE 

Having introduced the astrophysical sources and instrumen- 
tal effects we include in our simulations, we now summarise 
the method with which these are combined in order to pro- 
duce a full three-dimensional synthetic data set. This com- 
bination is performed in the Stokes representation and in 
image-space, however, some components, such as the cosmic 
signal, have been generated or modified or both in visibility- 
space (i.e., in the uv-plane of each frequency channel). The 
observation and instrumental parameters are initially spec- 
ified. These include the bandwidth, B, central frequency, 
vo, number of frequency channels, N c h, specific iro-tapering 
used and the parameters describing the Faraday screens. 

We begin by simulating and combining the astrophys- 
ical sources, starting with the cosmic signal, in terms of 
brightness temperature in a cube that corresponds to the 
specified observation. The EoR signal appears in Stokes I 
only and is simulated using the method described in Sec- 
tion [2T] The simulation snapshot redshift, z, and comoving 
side length, AD, of the EoR cube are given by 

z = ^ - 1 (64) 



and 



AD : 



c(l + zfB 
vixH{z) ' 



(65) 



where is the rest-frame frequency of the 21-cm line and 
H{z) is the Hubble constant calculated at redshift z. The 
EoR simulation is then instrumentally smoothed in the sky- 
plane for each frequency channel using the method described 
in Section l4~2l 

Next, the total intensity of DGSE is simulated as per 
Section |2.2.1| and its signal is added to Stokes /. Because 
we are assuming that there is no circularly polarised signal 
present, we set Stokes V to zero throughout. 



The first effect to be simulated is that of thin Fara- 
day screens. The corresponding polarised components of the 
DGSE are determined as per Section [2.2.2[ the plane of po- 
laristion is Faraday rotated and added to Stokes Q and U. 
For a single Faraday screen the intrinsic degree of polari- 



sation used in Equations ( 49 ) and ( 50 1 is simply the total 
brightness temperature of the polarised DGSE, i.e., fT^{y). 
This corresponds to the scenario where all of the polarised 
signal originates from behind the Faraday screen. In the sce- 
nario of multiple Faraday screens, the sum of intrinsic de- 
grees of polarisation for each screen equal the total bright- 
ness temperature of the polarised DGSE. Brcntjcns fc de| 



Bruyn ( 2005 1 provide an excellent illustration of more com- 



plicated emission and Faraday rotating scenarios. As previ- 
ously mentioned, for simplicity and to illustrate the princi- 
ple behind our cleaning method, we simulate the effect of a 
single thin Faraday screen only. 

Finally, the effect of instrumental polarisation leakage 
is simulated, using instrumental errors of the form discussed 
in Section 14.51 for each baseline. 



6 FOREGROUND SUBTRACTION 

In this section, we discuss the methods of removing diffuse 
Galactic foregrounds. By including only diffuse Galactic syn- 
chrotron emission we assume that all resolved discrete and 
extended sources have already been successfully removed 
from the observed maps, without any subtraction residuals. 
The requirements for subtracting point sources to a suffi- 



cient level are discussed by Liu et al. (20091, and a method 



for achieving this is presented by Pindor et al. ( Pindor et al 



2010) 



6.1 Continuum foreground subtraction 

There has been significant discussion in the literature re- 
garding low-frequency continuum foregrounds and strate- 
gies for their removal. While fluctuations due to these fore- 
grounds are expected to be several orders of magnitude 
larger than the reionisation signal (e.g., |Di Matteo et ak] 



2002 Oh & Mack 20031, the spectra of these foregrounds 



are anticipated to be smooth. As the reionisation signature 
includes fluctuations in frequency as well as angle, it has 
therefore been proposed that these foregrounds be removed 



through continuum subtraction (e.g., Gnedin & Shaver 2004 
Wang et al.[2006}[jehTet al.|2008| >, or using the differences 



in symmetry from power spectra analysis ( Morales & Hewitt 
|2004| |Zaldarriaga et al.|2004[ ). 

We model the foreground continuum using a polynomial 
of the form 



log 10 (T b G ) = 5>[log 10 0/)f 



(66) 



where v is the observed frequency and n is the highest or- 
der of log 10 (^) used for fitting. The appropriateness of the 
functional form of this fit can be illustrated by considering 
a pure DGSE signal along a single line of sight, with no in- 
strumental polarisation leakage or noise. In this case, we can 
perfectly fit our foreground model using a quadratic (n = 2) 



© 0000 RAS, MNRAS 000, 000-000 



12 Geil, Gaensler & Wyithe 



in \og 10 (u) with 

c o = logio(C) + a syn log 10 (Mo) - Aa syn [log 10 (!/ )] 2 , (67) 
Ci = -a syn + 2A« syn log 10 (Vo), 

C2 = -Aq s 



L syn : 



(68) 
(69) 

where C = f b G + <5T b G (f9). 

In the general case, we perform fits with this functional 
form along the line of sight for each spatial pixel in the 
simulated sky-plane. We then subtract the best fit, leaving 
residual fluctuations around the foreground emission with a 
near zero-mean brightness temperature, which will include 
any residual continuum foregrounds, leaked polarised signal, 
instrumental noise and the EoR signal. The consequence 
of removing the foreground continuum, together with the 
fact that interferometers do not make zero-spacing measure- 
ments, means that the true mean signal from each line of 
sight is unknown. As will be shown, this results in a de- 
creased contrast between ionised and non-ionised regions, 
and a loss of power from large-scale (small wavenumber) 
modes ( |Geil et al.|2008| >. 

6.2 Polarised foreground subtraction 



As discussed in Section |4.5| the leakage of polarised fore- 
ground components into the Stokes I signal can lead to sig- 
nificant contamination of the cosmic signal. The procedure 
for removing this contamination is shown schematically in 
Figure [6] For each line of sight: 

(i) the observed complex linear polarisation, P(A 2 ), is 
Fourier transformed to find its corresponding Faraday dis- 
persion function, F(<f>); 

(ii) the number of peaks, Nfs, is determined from F(<f)) 
by a thresholding technique or a peak removal procedure 
similar to the method discussed here (F is used rather than 
Fj since the signal-to-noise ratio is greater and the polarised 
signal is 'pure' in F), and each peak's Faraday depth and 
intrinsic polarisation angle is found using Equation Q; 

(iii) in the case that Nfs ~ 0, there is no polarised con- 
tamination to clean; otherwise, 

(iv) we attempt to remove each identified peak from the 
Stokes / signal. This is done by Fourier transforming /(A 2 ) 
to give the effective Stokes I Faraday dispersion function, 
Fi((f>), then removing the largest peak by constructing a 
model of the real and imaginary parts of a polarised compo- 
nent (using the <j> an d Xo found from the analysis of F and 
leakage solutions found by calibration); 

(v) the cleaned Fi{cf>) is inverse Fourier transformed back 
to give the cleaned Stokes I signal. 

This process is repeated for each line of sight, and the 
cleaned Stokes / signals are reconfigured into a three- 
dimensional data cube. It is possible to remove polarised 
foreground components in a number of iterations using a 
loop-gain factor 7 to bring the peaks down slowly, treat- 
ing the strongest component in each iteration. This method 
is similar in spirit to the CLEAN deconvolution algorithm 
| Hogbom|1974 l used in radio imaging. However, we subtract 
a peak in a single iteration in our simulations. 

The single iteration procedure described above performs 
well for lines of sight where the number of Faraday screens is 
correctly calculated. This is the case when the RM peaks are 




/ - 




> 




Remove largest 
peak from Fi 




Figure 6. Flow diagram for single line-of-sight polarised fore- 
ground removal. 



spaced sufficiently far apart in <^)-space such that one does 
not superimpose another (the differentiation between peaks 
also depends on the rotation measure spread function, R(4>), 
of the instrument at a particular frequency). The single it- 
eration procedure fails when the different Faraday screens 
are close to each other, however, follow-up analysis of lines 
of sight, including higher dimensional fitting, where this is 
an issue can be performed. 



6.3 Summary of assumptions and other 
considerations 

The aim of this paper is to present the important first 
step toward solving the important problem of polarised fore- 
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ground leakage. In order to do this a number of assumptions 
have been made throughout this work so as to simplify the 
simulations and attempt to demonstrate the basic principal 
of our algorithm. 

(i) We have only simulated a single thin Faraday screen as 
a means to simply demonstrate our idea of leakage removal. 
The real situation is of course more complicated, including 
multiple thin screens and/or slabs, but we hope to tackle a 
more realistic scenario in later work. 

(ii) The effect of the instrumental beam has been ne- 
glected. Any effect of the beam should of course be well 
understood and modelled. We assume that these effects have 
been treated prior to the application of our procedure. 

(iii) There are many contaminants and foregrounds to 
consider, however, any components not mentioned already 
have been considered insignificant to the present procedure. 

(iv) We have not made a realisation of a power-spectrum 
for the spatial fluctuations in non-polarised foreground. A 
Gaussian field has been used for the sake of simplicity. 

(v) Observations in our simulations are assumed to have 
been made in the zenith direction. In a more realistic model 
observations away from the zenith would need to be consid- 
ered and the impact of dipole projection considered. 

(vi) The size of the errors we consider is no doubt larger 
than what be actually be the case. This was chosen in order 
to make a clear demonstration of the leakage. Also, we have 
modelled the errors with no bias despite the possibility that 
in reality a systematic bias may be present due to the mass- 
manufacturing method of the dipoles. 

(vii) We make an important assumption that Faraday 
rotation effects of the ionosphere have already been ac- 
counted for. We argue that these ionosphere through which 
the instrument is observing is well measured and understood 
through the use of polarised point-source callibrators and 
that any effect has been removed (by our method or other- 
wise) prior to the application of our procedure on Galactic 
Faraday rotating regions. 

Naturally, at some point when the level of noise is simi- 
lar to the level of polarised foreground leakage in the Faraday 
dispersion function a confusion limit will be found. This has 
been tackled in the following way. First, by fitting for RM 
and polarisation angle using F (rather than Fi) where the 
level of confusion will be less since the peaks correspond- 
ing to Faraday-rotating structures will be larger than their 
Stokes I counterparts. Secondly, since the data cubes com- 
prise of integrated data, the level of noise can, to a certain 
degree, be controlled by the integration time. 



7 SPHERICALLY AVERAGED POWER 
SPECTRA 

In this paper, we employ the dimensionless power spectrum, 
A§i(fc,2), as the key statistical measure, which is the con- 
tribution to the variance of the redshifted 2f-cm bright- 
ness temperature contrast, STh(z), per logarithmic interval 
in wavenumber, k = 2n/X. This measure is related to the 
dimensional form of the power spectrum, P2i(fc, z), by 



P2i(fc, z) is estimated by averaging over all m modes of the 
Fourier transform (T) of T(z) in a thin spherical shell in 
fc-space, 



1 m 

P 21 (k,z) = {\f(k,z)\ 2 ) k = i£)|f 4 (fc,*)| a 



(71) 



Note that we must use caution when employing this statis- 
tical measure since the spherical symmetry of the signal is 
broken by redshift evolution over certain ranges of redshift 
(Az w 2 for a 32 MHz band centred on z = 8.5). This is- 
sue has been discussed by a number of authors (see, e.g., 
Morales fc Hewitt ||2004l |Barkana fc Loeb|[2005l |McQuinn| 
et al.| 2006 I and its effect on the sensitivity of the 21-cm 
power spectrum has been investigated by [McQuinn et al.| 
l |2006[ ). Following the removal of both the continuum and po- 
larised foregrounds, the full 32 MHz data cube is divided into 
four 8 MHz sub-bands (Az <; 0.5 for z ov < z < 8.5), each of 
which is independently reduced into a power spectrum. This 
sub-banding, together with continuum foreground removal, 
eliminates measurements of the largest spatial scales along 
the line-of-sight direction (these corresponding to the lowest 
fc|| modes) and imposes a minimum accessible wavenumber 
of fc min w 0.04[(1 + ZV7.5]" 1 Mpc- 1 ]. 



7.1 Sensitivity to the 21-cm power spectrum 

Our numerical power spectrum measurements are subject 
to sample variance, arising from the finite number of inde- 
pendent modes counted in each fc-shell, which corresponds 
to the finite number of independent wavelengths that can 
fit into the simulated volume. Any radio interferometer is 
also subject to instrumental noise and has limited sensitivity 
arising from the finite volume of the observation. We con- 
sider observational parameters corresponding to the design 
specications of the MWA, and of a hypothetical follow-up 
to the MWA (which we term the MWA5000). The MWA 
5000 is assumed to follow the basic design of the MWA. The 
quantitative differences are that we assume the MWA5000 
to have 5000 tiles within a diameter of 2 km, with a flat an- 
tenna density core of radius 40 m. In each case, we assume 
one field is observed for an integrated time of 1000 hr. 

We compute the sensitivity with which the 21-cm power 
spectum could be detected following the procedure outlined 



by McQuinn et al. (20061 and Bowman et al. (20061. Sec- 



tion |A2| of the Appendix contains a full derivation of the 
sensitivity-related quantities used. Written in terms of the 
cosmic wavevector fe = kit + k±, where feii and k± are the 
components of k parallel and perpendicular to the line of 
sight respectively (see Appendix |A2[ ), the resulting error in 
the 21-cm power spectrum per mode is 



5P 2 i(k) 



r s 2 y8 D 2 AD 
Bt n h {U,v) 



+ P 2 i(fe), 



(72) 



A|r(M) 



(2n) 



-4irk 3 P 2 i(k,z) 



(70) 



where T ByB ~ 250[(1 + z)/7] 2 ' 6 K is the system tempera- 
ture of the instrument, D(z) is the comoving distance to the 
point of emission at redshift z, AD is the comoving depth 
of the survey volume corresponding to the bandwidth B, to 
is the total integration time, nh(U,v) is the number den- 
sity of baselines that can observe the visibility U, where 
U = k±D/2ir and A is the observed wavelength. 

Although the observed 21-cm power spectrum is not 
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spherically symmetric, it is symmetric about the line of 
sight. This makes it possible to calculate the overall power- 
spectral sensitivity of the radio interferometer using the 
Fourier modes contained within an infinitesimal annulus 
around the line of sight of constant (k,6), where cos(#) = 
k ■ z/k [z is the unit vector pointing in the direction of the 
line of sight). The power spectral sensitivity over such an 
annulus is given by 



a P (k,8) 



6P2i{k,6) 
y/N m {k,e)' 



(73) 



where N nl (k, 8) denotes the number of observable modes in 
the annulus (only modes whose line-of-sight components fit 
within the observed bandpass are included). In terms of the 
fc-vector components k and 8, the number of independent 
Fourier modes within an annulus of radial width dk and 
angular width d8 is N m (k,0) = 2nk 2 V sin(8)dkd8/(2n) 3 , 
where V = D 2 AD(X 2 /A c ) is the observed volume. Averag- 
ing <jp(k,8) over 8 gives the spherically averaged sensitivity 
to the 21-cm power spectrum ap(k), 



,(*) 



E 



>(M)' 



(74) 



In order to find the sensitivity in terms of Afi we use Equa- 



tion (70 1, which gives 
1 



a A 2(k, z) = 



(27T) 



-47T/C 



,{k,e,z) 



-1/2 



(75) 



of these components {dashed) are shown for comparison. 
The third row of plots show the cleaned effective Stokes / 
Faraday dispersion function with the contaminated effective 
Stokes / Faraday dispersion function [faint) shown for com- 
parison. The bottom row of plots show the cleaned signal 
(black) after being Fourier inverted and the model cosmic 
signal (blue) for reference. 

It is evident that the instrumental polaristion leakage 
due to the Faraday screen at the smaller depth (shallow 
Faraday screen) is less effectively removed than the deeper 
screen. This is due in part to a residual long- 'wavelength' 
fluctuation in the signal after continuum foreground removal 
which leaves small-0 power in the effective Stokes I Fara- 
day dispersion function. Residuals such as these can make 
a small-<^ peak in the Stokes / Faraday dispersion function 
asymmetric and therefore more difficult to fit accurately. 
Deeper peaks are not affected in this way as the modelled 
continuum foreground does not contain fluctuations that 
masquerade at these Faraday depths, therefore no residu- 
als at these depths are present after continuum subtraction 
making more accurate fitting possible. 

Unwanted small-<j!> residuals have been removed by ap- 
plying a high-pass filter to the cleaned effective Stokes / 
Faraday dispersion function. The extent of this filter is 
^ (f>{ <. 21.2 radm -2 , which has an equivalent width in 
frequency-space of approximately 4.7 MHz for the correlator 
specifications used here. Negligible cosmic signal is removed 
by this procedure since the EoR fluctuations are of shorter 
'wavelength' at the corresponding redshift at this frequency. 



8 RESULTS 

In this section, we present the results of our simulations 
and analysis. Section |8.1| shows the result of applying the 
polarised foreground removal algorithm for a single line of 
sight. Section [8. 2| gives the results for a full three-dimension 
data set. Figure [7] shows the resulting spherically averaged 
21-cm power spectra. 

8.1 Single line of sight 

Figure[7]demonstrates the removal of a single polarised fore- 
ground component from a noiseless Stokes / signal. For 
this example, we assume that any continuum DGSE has 
already been removed, leaving only the cosmic signal and 
leaked polarised foreground. Two different values of Fara- 
day depth have been modelled in order to show the differ- 
ence in cleaning performance. The left-hand example uses 
</> = 35.5 radm" 2 with p(A§) w 176.8 K and xo » 0.794. 
The right-hand example uses 4> = 5.5 radm -2 with p(Aq) « 
161.6 K and xo ~ -1.52. 

The top Stokes / visibility versus frequency plots in Fig- 
ure [7] show the signal to be cleaned (black) and the model 
cosmic signal (blue) for reference. Note that the frequency 
axis is linear and therefore the apparent 'period' of the con- 
taminated signal containing the leaked polarised foreground 
is non-uniform. This is because of the A 2 -dependence of the 
plane of polarisation as per Equation jsj. The second row 
of plots show the effective Stokes / Faraday dispersion func- 
tion, Fj((j>), of the contaminated signal, including the real 
and imaginary part as well as its modulus. The fit to each 



8.2 Full data cube 

Figure [8] shows a simulated three-dimensional Stokes I 
data cube over a 32 MHz band centred on 157.8 MHz 
(6.3 <C z <. 7.8) using 128 channels. This simulation in- 
cludes both continuum and polarised foreground removal. 
Faraday depth has been modelled as a gradient over sky- 
plane with 1 ^ (f> ^ 5 radm -2 and the intrinsic polarisa- 
tion angle for each line of sight has been randomly assigned, 
— 7r/2 < x ^ 7r /2 rad. The upper-left cube shows the full 
Stokes I signal before cleaning. The upper-right cube has 
had the continuum DGSE fitted and removed in each line 
of sight, leaving the cosmic signal and the effects of instru- 
mental polarisation leakage. The banding along the 8 X - and 
^-frequency planes result from the oscillatory behaviour of 
the leaked polarised signal. The middle-left cube is the raw 
fully cleaned signal having had the leaked polarised signal 
removed. The pre-leakage-cleaned banding has mostly dis- 
appeared, with the exception of channels close to the edge of 
the frequency band. This is a systematic effect due in part to 
the continuum foreground removal. The central-right cube 
shows the fully cleaned signal after high-pass filtering of the 
effective Stokes / Faraday dispersion function of each line of 
sight. The same filter width has been used as for the one- 
dimensional example in Section |8.1| The central frequency 
slice of the fully cleaned and high-pass filtered cube is shown 
at the lower-right, with the model EoR signal shown at the 
lower-left for comparison. The model EoR cube has not had 
each line of sight's mean subtracted, which would occur as 
a consequence of the continuum foreground subtraction. 

Figure [9] shows a selection of line-of-sight signal profiles 
from the three-dimensional simulation for the model cosmic 
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Figure 7. Examples of the removal of a single polarised foreground component from Stokes /. Top: Synthetic Stokes / signal (black) 
comprising of EoR signal (blue) and contamination due to instrumental polarisation leakage. Second from top: Effective Stokes / Faraday 
dispersion function (solid) of total signal and fits to its magnitude and real and imaginary parts, (dashed). Second from bottom: Residual 
effective Stokes I Faraday dispersion function after cleaning. Original effective Stokes I Faraday dispersion function (faint) is shown 
for comparison. Bottom: Cleaned Stokes I signal (black) and EoR signal with line-of-sight mean removed (blue) for comparison. Left: 
Single Faraday screen at <f> = 35.5 radm -2 with p(Aq) Ri 176.8 K and xo ~ 0.794. Right: Single Faraday screen at <j> = 5.5 radm -2 with 
P(Aq) 161.6 K and xo ~ —1-52. A high-pass filter has been applied to the cf> = 5.5 radm -2 example in order to eliminate the residual 
long-'wavelength' fluctuation that remains after cleaning (faint dashed). The filter edge at <f>{ fa 21.2 radm -2 is shown by the vertical 
line on the plot of the cleaned effective Stokes / Faraday dispersion function. 



signal (black), fully cleaned signal with no high-pass filtering 
of their effective Stokes / Faraday dispersion function (red) 
and fully cleaned signal after high-pass filtering of their ef- 
fective Stokes / Faraday dispersion function (blue). 



8.3 Power spectrum measurement in the presence 
of foregrounds 

Here, we summarise the effects of various foreground clean- 
ing methods on the potential power spectrum measurement. 
Figure [To] shows the spherically averaged 21-cm power spec- 
tra for a foreground cleaning run performed on a simulated, 
noiseless Stokes I signal over a 32 MHz band, centred on 
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Figure 8. Simulated noiseless Stokes / signal over a 32 MHz band, centred on 177.5 MHz (6.3 < z <^ 7.8) using 128 channels. Faraday 
depth has been modelled as a gradient over sky-plane with 1 ^ <f> ^ 5 radm - 2 . Intrinsic polarisation angle for each line of sight has been 
randomly assigned. Top-left: Full noiseless Stokes I signal showing domination by DGSE. Top-right: Stokes I signal after continuum 
foreground removal. Middle-left: Stokes / signal after continuum and polarised foreground removal. Middle-right: Stokes / signal after 
continuum and polarised foreground removal and high-pass filtering. Bottom-left: Model EoR signal. Bottom-right: Central channel slice 
through Stokes / signal after continuum and polarised foreground removal and high-pass filtering. 
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175 180 
v [MHz] 

Figure 9. Simulated noiseless Stokes / signal over a 32 MHz band, centred on 177.5 MHz (6.3 < z <C 7.8) using 128 channels. Faraday 
depth has been modelled as a gradient over sky-plane with 1 ^ ^ 5 radm -2 . Intrinsic polarisation angle for each line of sight has been 
randomly assigned. Example linc-of-sight signal profiles: model EoR signal with line-of-sight means removed (black), fully cleaned signal 
with no high-pass filtering of their effective Stokes I Faraday dispersion function (red) and fully cleaned signal after high-pass filtering 
of their effective Stokes I Faraday dispersion function (blue). 



177.5 MHz (6.3 <; z <; 7.8) using 128 channels. Instrumen- 
tal sensitivity curves for MWA500 (dashed) and MWA5000 
(dot-dashed) and 1000 hr integrations are also shown, where 
it has been assumed that the total 32 MHz bandwidth has 
been split into four 8 MHz sub-bands. 

It can be seen that the power spectrum of the cleaned 
signal agrees very well with the power spectrum of the model 
cosmic signal, once it has been apodised and further cleaned 
with a high-pass filter, within the sensitivity limits of the 
MWA. The difference between the two for k <;0.04 Mpc -1 
is due to the loss of large-scale contrast owing to continuum 
foreground subtraction. The up-turn at large k (small scales) 
arises from the residual fluctuations present at the edges of 
the bandpass (as a result of polarised foreground removal). 

Apodisation of the bandwidth is required in order to 



eliminate spurious power in the power spectra due to the 
band-edge effects discussed in Section [8. 2| No overall power 
is lost through this procedure, however, because the sampled 
volume is smaller, there will be a decreased effective sensi- 
tivity due to decreased number of Fourier modes, 7Y m oc V 
[see Equation ( 73 1] . 



9 SUMMARY 

In this paper, we have shown how polarised signals can 
contaminate the non-polarised signal — and therefore con- 
taminate the signature of cosmic reionisation — through the 
process of instrumental polarisation leakage. We have also 
demonstrated how, in principle, RM synthesis may be used 
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Figure 10. Spherically averaged 21-cm power spectra. Beam-smoothed model EoR signal; Stokes / signal after continuum and polarised 
foreground removal, high-pass filtering and band apodisation at redishifts 2 = 7. Instrumental sensitivity curves for MWA500 (dashed) 
and MWA5000 (dot-dashed) and 1000 hr integrations are also shown, where it has been assumed that the total 32 MHz bandwidth has 
been split into four 8 MHz sub-bands. 



to recover the cosmic signal by applying an RM synthe- 
sis technique to a three-dimensional synthetic data cube. 
In demonstrating this, we have assumed that the brightest 
point sources have been sufficiently removed, and modelled 
the brightest astrophysical sources of foreground contami- 
nation only, for the simplest case of a single, thin Faraday 
screen. 

Our results show that it is possible to accurately recover 
the signature of reionisation late in the epoch of reionisation 
(z ~ 7), by way of the 21-cm power spectrum, within the 
sensitivity limits of an instrument similar to the MWA. In 
addition to the statistical signature of cosmic reionisation, 
we have also shown that the tomographic image of the per- 
colating IGM can be successfully cleaned of contamination. 

We have applied our technique using a simple fore- 
ground model and a single thin Faraday screen. We believe 
it is possible to extend the application of this technique to 
more complex and realistic foreground and Faraday screen 
models by performing the cleaning process iteratively. This 
is an important improvement for lines of sight that have 
multiple Faraday-rotating screens or screens that are thick 
enough to be resolved by the instrument. For situations 
where a line of sight contains multiple screens, well sepa- 
rated in 0-space, the single iteration cleaning method used 
in this paper may be applied to each foreground peak in 
the effective Stokes I Faraday dispersion function. Although 
contamination of the Stokes / signal due to very shallow 
Faraday screens (which have fluctuations in frequency-space 
on a scale larger than the EoR signal) can be removed by 
simply applying a high-pass filter in 0-space with minimal 
effect on the recovered cosmic signal, Faraday screens that 
lead to contamination that fluctuate on scales similar to the 



EoR signal (large \4>\) cannot, but can be easily removed 
using our method. 

It has already been shown that neither point sources, 
continuum DGSE nor instrumental noise will pose insur- 
mountable obstacles to the detection of the cosmic 21-cm 
signal (see e.g., | Wang et al.|[2006] |Geil et al.||2008| [Jeli6| 



et al.|2 008 Bo wman et al.|2009||Liu et al.|2009^ Our results 
suggest that the cleaning of leaked polarised foregrounds will 
also be possible. Removing this contaminant will be a critical 
element in the removal of foregrounds that will be essential 
in order to detect the cosmological 21-cm signal. 
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APPENDIX A: STATISTICAL ERRORS AND 
INSTRUMENTAL SENSITIVITY 

Al Fiducial array configuration 

A 1.1 Coordinates and notation 

The visibility measurements of the next generation of low- 
frequency radio arrays will sample the instrumental response 
from an inherently three-dimensional volume of space at 
high redshift. Only the fluctuating part of the signal con- 
tributes to an interferometric visibility, which we denote 
by AV(u,v, v), where u and v have their usual meaning 
(see, e.g., Thompson et al."||2001 |. As discussed by |Morales| 
( |2005[ ), there are three useful representations of the data: 
the image representation, the visibility representation, and 
the full Fourier representation. The image cube is coordi- 
nated by sky position and frequency (9 x ,0 y ,u) and can be 
transformed into a cube of visibilities by performing a two- 
dimensional Fourier transform in each frequency channel. 
Instrumental visibilities have coordinates given by (u,v, v), 
where u = 1/8 X and v = l/9 y . The visibility cube is trans- 
formed into the full Fourier representation by performing a 
one-dimensional Fourier transform in frequency, giving coor- 
dinates in terms of (u, v,rj), where r\ = \ jv is the wavenum- 
ber corresponding to the frequency of the visibility mea- 
surement. Here, we preserve the traditional meaning of the 
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interferometric quantities u and v by defining the three- 
dimensional wavevector of a Fourier mode by 



T = ux + vy + rjz — U + r/z. 



(Al) 



The cosmological Fourier mode k = k± + fen and its rela- 
tionship to the instrumental Fourier mode T is discussed in 
Section |A2l 



A1.2 Antenna distribution - continuous approximation 

Although the number density of antennas for an interferom- 
eter will be discrete in practice, here we make the simplifying 
assumption of a continuous distribution. The continuous an- 
tenna distribution generalises our analysis (to some limited 
degree) and is justified for large numbers of antennas in a rel- 
atively compact configuration (i.e., a sensible distribution, 
considering the operating frequency, thereby providing full 
wv-coverage). For a circularly symmetric array of AT a anten- 
nas (tiles) with a constant antenna density core of radius r c , 
maximum radius r max and an inverse square antenna density 
profile for r > r c , the properly normalised antenna number 
density as a function of r is given by 



r <r c 



n a (r) 







— r c < r < 

' max 



(A2) 



where n c is the antenna number density of the core, given 
by 

Na. 



7rr c 2 [l + 21n(r max /r c )]' 



(A3) 



This is capped such that there is a maximum density of 
one antenna per 36 m 2 . The baseline number density can be 
written 



2nrdr n a (r) / #n a (|r - XU\), (A4) 



n h (U,v) = C h 



where C = C\>{y) is a frequency dependent normalisation 
factor such that / dUn h {U,v) = N^N* - l)/2. The fre- 
quency dependence of rib arises due to the frequency depen- 
dence of the baselines (U otv). 

Figure |A1| shows the corresponding antenna number 
density for a continuous n a oc r~ 2 distribution, and the 
continuous baseline number density for v = 80, 158 and 
300 MHz. 



A2 Power spectrum measurement 

Calculations of the sensitivity to the 21-cm power spectrum 
for an interferometer have been presented by a number of 



authors ( 


Bharadwaj & Sethi 


2001 


|Bharadwaj & Pandey 


2003 


Morales 2005; |Bowman et a 


.(|2006| |McQuinn et al. 


2006 


1. W 


3 follow the procedur 
•awing on results fron 


e outlined by 


McQuinn et al. 


(20061, di 


1 Morales ( 2005 1 and Bowman 



et al. (2006) for the dependence of the array antenna density 
on radius, n a (r). 

As mentioned in Section [A 1| the baseline number den- 
sity rib is a function of frequency. This dependency results 
in a smaller range of baseline values and denser sampling 
at lower frequencies compared with the coverage at higher 
frequencies (assuming an identical visibility regridding scale 
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Figure Al. Top: Antenna number density for a continuous n a oc 
r~ 2 distribution. Bottom: Continuous baseline number density for 
v = 80, 158 and 300 MHz, corresponding toz= 16.8, 8 and 3.7 
respectively. 



is used at each frequency). We ignore this dependency here 
since for bandwidths of B <g 32 MHz, centred upon a fre- 
quency of interest with respect to the EOR, the relative scale 
difference in baseline for any antenna pair at each end of the 
frequency band is 0.8 <; U (v m i n ) /U (zwx) = l-B/v ma . K < 1. 

The uncertainty in a measurement of the power spec- 
trum has two separate components, due to (i) the thermal 
noise of the instrument and (ii) sample variance. Here we 
introduce the power spectrum in the full Fourier represen- 
tation and the first-order covariance matrices used to de- 
rive the relevant measurement errors, as well as discuss the 
relationship between instrumental and cosmological wave 
modes. Instrumental noise and sample variance are discussed 
in Section |A2.1| and Section |A2 . 2| respectively. 

The wavevector k can be conveniently split into two 
orthogonal components, one parallel to the line of sight and 
the other transverse to it, such that fe = feii + k±. The 
transverse wavenumber k± — \k±\ is given by 

c dz' 



2tt 



2nU 



D{z) 



(A5) 



k = — = 

± d ± D(z)A6 D(z)' v ~' Jo H(z' 

where d±_ is the comoving dimension of the mode perpen- 
dicular to the line of sight, A9 is the angle subtended by d± 
on the sky (in radian) and D(z) is the comoving distance to 
the point of emission at redshift z. The parallel wavenumber 
fen = |feii| is given by 



2tt 



2ttB 
ADAv 



AD 1 



c(l + z) 2 B 

V2lH{z) ' 



(A6) 



where dii is the comoving dimension of the mode parallel to 
the line of sight, AD is the comoving depth of survey volume 
corresponding to the bandwidth B and 1/21 is the rest-frame 
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frequency of 21-cm line. dU = [/ max /iVjy and dr\ = 1/B with 

t)max = N ch /B. 



dk± = 



2-wdU 



dkw = 



2tt 
AD' 



(A7) 



Using these definitions we can write P2i(fc) = 
(D 2 AD/P)P 2 i(Y). 

In the full (u, v, rf) Fourier representation, the measured 
visibility V(U , v) is given by 



J(T) = / dvV(U,v)exp{2nivr)), 



(A8) 



where the tilde indicates the two-dimensional Fourier trans- 
form that relates the specific intensity I(0 X , 8 y , v) in the im- 
age cube with the visibility V(U, v) in the visibility cube. 

The error in the 21-cm power spectrum per mode of 
volume d 3 T = dudvd-q fti A C /X 2 B is given by 

<5P 2 i(Y) = [c N (T) + C sv (T)]d 3 T, (A9) 

where C N (Y) is the rms of instrumental noise and C SV (T) 
is the rms of the sample variance. The covariance of two 
measured I(T) values is 

Cf 3 = C x (Ti, Tj) 



(A10) 



where C X (T) = C X (Y, Y) and X represents the source of 
the error. 



instrumental noise is found using Equation (AlOl, and is 
I (T 1 )[/ N (Y J )] 



a 



^ V^ na {U i ,v m )exp(2niv m r li )A 

m = l 

V r r m S (fi, f n ) exp(27rii> , n77j)A 

,n=l 

JV ch 
m,n — 1 

= (A;/) 2 £ (v£,(Ut, u m ) [v^U^Vm)] 



m — 1 
m — l 

/■'N ITT ,,\ — T/N 



(A14) 



where Vj^^U, v) = V TU1B (U, v) exp(2-Kivr\) . In the case of 
a contiguous bandwidth with B <C fo, KmsC^i,") is ap- 
proximately the same in all frequency channels and the last 



expression in Equation ( A14 1 simplifies to 
C% = N^AvflvrUU^vrfS^ 



a 2 t s , 

\ 2 BT B , 



B 

2 



BA c t n h (U t ,Vi) l3 ' 



(A15) 



A2.1 Instrumental noise 

The rms noise fluctuation per visibility per frequency chan- 
nel is given by ( jRohlfs fc Wilson|2004| | 

9T b 2fc B T sys 



oI v A C \J Avtjj 



K. 



(All) 



where T sys is the system temperature (K), A e is the effective 
area of one antenna (m 2 ), tjj is the integration time for that 
visibility (s), Av is the frequency bin width (Hz) and fca 
is the Boltzmann constant. For the system temperature, we 
use the mean sky brightness due to DGSE in Equation pi). 
Instrumental noise is uncorrelated in both U and v, i.e., 
{V^(Ui,v m )\Vi£ui(Uj,Un)]*) oc Sij5 mn . The average inte- 
gration time tu that an array observes the visibility U is 
given by 



tu = -^t n h (U, u), 



(A12) 



where to is the total integration time and nh(U, v) is the 
number density of baselines that can observe the visibility 
U. 

Performi ng t he discrete version of the transformation 
in Equation ( A8 1 for the noise-only rms visibilities ^ in 



Equation (All l gives 



7 N (T) = ^ V^[U,v m )exv(2nw m rj)Ai 



(A13) 



where iVch = Bj Av is the number of frequency channels, 
each centred on a frequency Vj. The covariance matrix for 



A2.2 Sample variance 

The second component of uncertainty is due to sample vari- 
ance within the finite volume of the survey. The covariance 
matrix for sample variance is given by 



= I 



21 (Yi) J /2 l (Yj) j* 



d J YP2i(Y)w(Y l - r)w*(r, - T) 



\ 2 B 



P2i(Ti)<5 i: 



(A16) 



where W is the instrumental window function given by 
the field of view and bandwidth of the array. The last line 
follows by performing the integral after noting that the 
Fourier transform of the window function, W is non-zero 
only in a volume dudvdn « A C /X 2 B and therefore W is 
normalised by (dudv drj)' 1 ^ 2 . 



Combining Equations ( |A9| , ( |A15[ ) and ( |A16[ ), the re- 
sulting error in the 21-cm power spectrum per mode is 



<5P 2 l(Y) — / -*sys 



1 



+ P 2 i(T). 



A c J torih{U,v) 
Written in terms of the cosmic wavevector, 



AP lu\ P S ye D 2 AD f\ 2 ^ 2 
where U = k±D/2n. 



+ P 2 i(fe), 



(A17) 



(A18) 
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A3 Uniform gains 

For a single baseline with uniform gains (equal to unity, 
therefore all 7's vanish), Equations (|61[)-( [63"| ) simplify to 

I'ii = I ij + \8 + -U ij -\i8-+Vi h (A19) 

Q'ii = Qa + - \i5-Va, (A20) 

Vii = U lJ + U + .I lJ -h ++ Q l] . (A21) 

Mij has to be summed over all baselines and their conju- 
gates however, which then gives the following expressions for 
the response in Stokes /, Q and U (all of which are frequency 
and sky position dependent): 



(A22) 



\ Vb,6* V6,b* 

U' = U+^(lJ25+--Ql2s+\ (A24) 

V V6,b* V6,b* / 

where A<b = N(N — l)/2 is the number of baselines and 
I = y^I tj /(27Vb) etc. If we now assume that the sources have 
no circularly polarised component (V ~ 0) and ^ 5R[A^,,] ~ 
~}2 $l[Ayy] (this constraint only affects the Stokes I leakage 
into Stokes U), where the A i terms are defined by 

^xx — di,xy i dj xy, (A25) 

Aj.^ = di tXy zb dj yx , (A26) 
A ya; = di^yx ^dj X y, (A27) 
A M!/ = di t y X i dj yx , (A28) 



then Equations (A22|-(A24l simplify to 
T = /, (A29) 

Q' = <5+2^X)X)[»(A+)+»(A+)], (A30) 

u' = ^^EE^ + ^i]- ( A31 ) 

The linearly polarised components Q and U can be written 

2 Mi L r)(g)> 

where Lq^ and Luq are real, frequency-dependent Stokes 
leakage parameters. Equation ( |A32[ ) has been used in the 
example in Figure [5] which shows simulated examples of 
the CHECK Faraday dispersion function, F((j>), for a single 
baseline illustrating the effect of leakage between Stokes Q 
and U, noise and DC offsets. In general, however, the 
gain/leakage errors for each baseline will be independent, 
thereby complicating the analysis above. 
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